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SUMMARY 

An equation of motion about a satellite spin axis is derived 
that includes solar radiation torques produced by a solar paddle 
array. The equation is applied to the Explorer XII and solved to 
explain the observed increase in the satellite spin. Also, an ap- 
proximate equation is developed which permits rapid calculation 
of the effect of solar torques on satellite spin. 
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THE EFFECT OF SOLAR RADIATION PRESSURE 
ON THE SPIN OF EXPLORER XII 


by 

J. V. Fedor 

Goddard Space Flight Center 


INTRODUCTION 

Recent Goddard Space Flight Center satellites have had relatively large solar cell paddles to 
obtain power from the sum An interesting effect of this has been the marked change of satellite spin. 
For the Explorer XII (launched August 16, 1961), the spin increased approximately 20 percent from 
its initial value in four months. It is the purpose of this note to report an analysis which shows that 
the increase of the Explorer XII spin can be attributed to solar radiation pressure. 


ANALYSIS 

It is a well established fact that sunlight is electro-magnetic energy. From Einstein*s equa- 
tion E = me 2 , where E is energy, m is mass, and c is the velocity of light, a mass and thus a 
momentum can be associated with sunlight. Whenever light impinges on a surface (for example, 
solar paddles on a satellite) there is a change of momentum and hence a force (torque) developed 
which in the case of spinning satellite may effect the rotation. Although this torque is usually small, 
it can under certain conditions produce significant results. 

In the interest of reducing algebraic complexity, improving clarity, and to develop equations for 
engineering application, a restricted point of view will be taken in the subsequent analysis. That is, 
motion and torques only about the satellite spin axis will be considered. For the Explorer XU, this 
is found to be adequate to describe spin changes. A more general and abstract approach to solar 
torques acting on a satellite is given in Reference 1 and 2. 

Once the satellite is injected into orbit and any nutation angle has been damped out due to energy 
dissipation, the satellite can essentially be considered as a one degree of freedom system (A favor- 
able moment of inertia ratio is assumed). The effect of forces other than around the spin axis usu- 
ally can be neglected if satellite spin is relatively high (20 to 30 rpm with a spin moment of inertia of 
about 3 to 4 slug feet 2 ). Moments tending to affect the spin axis are mainly magnetic, solar radia- 
tion, and aerodynamic. Because of the highly eccentric orbit of the Explorer XII and relatively high 
perigee (after one month in orbit: apogee 77, 199 km, perigee 448 km), the aerodynamic moments 
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are small and can be neglected. The Explorer XII was also designed to minimize magnetic torques 
and the high apogee tends to reduce magnetic effects even further. Hence, magnetic torques will also 
be neglected. Thus, only solar radiation torques are left to be considered. 

Figure 1 shows the general configuration of the satellite. Only direct radiation torques produced 
by the solar paddles will be considered; radiation torques produced by the main satellite body geometry 
about the spin axis are considered negligible. The equation of motion around the spin axis is simply 

^TORQUES = I fjr (1) 

where o - dO/dt, 0 is the angular coordinate around the spin axis, and I is the spin axis moment of 
inertia. 

In developing the moment expression it will be assumed that there are two distinct pairs of pad- 
dles. The force acting on a paddle is the product of the radiation pressure factor, P(lbs/ft 2 ), and 
the projected area perpendicular to the sun T s rays, A(ft 2 ): 


Incident light force - PA. 


For the reflected light, a reflectivity coefficient, p, must be included 

Reflected light force “ pP A . 

In general P varies with the earth-sun distance and this is discussed later. The reflected light from 
a solar paddle is a rather complicated phenomenon being dependent upon the optical properties of 
the reflecting surface and the impinging light. It is assumed here that the reflection is spectral, that 
is, angle of incidence equals angle of reflection and that the reflectivity is independent of the incident 
angle. Figure 2 shows the geometry for the torque equation. The partial component of torque about 
the spin axis for the incident light is 


PAh | cos 0 | sin 0 , 


where h is the distance from the spin axis to the centroid of the paddle, and <$> is the sunline-spin 
axis angle. For the reflected light it is 


pPA | h cos 8 j sin (</> - 2ns) 


where ns is the angle between the normal to the paddle and the direction of the sun. Absolute signs 
are used on the cosine term since the moment arm is always positive. Summing the torques for each 
paddle and considering counterclockwise positive, we obtain 
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£ TORQUES = - P A ! h i I cos 5 | (sin <p + p sin [cp - 2ns,)j 

+ PA 2 h 3 | sin 0 | (sin </> + p sin (<£ ~ 2ns 2 )^ 

+ PA 3 hj | cos 0 | (sin <p + p sin (<£ - 2ns 3 )j 

- PA 4 h 2 | sin 0 | (sin <p + p sin [<p - 2ns 4 )^ 


( 2 ) 


Equation 2 is the basic torque expression for direct solar radiation effects about the spin axis due to 
the solar paddles. 

The projected area of a paddle as a function of <f>, 8, and the geometry of the paddle is developed 
in Appendix A. Using the basic projected area equation, the projected area of the four paddles can 
be expressed in the following way 


Z-AXIS 




Figure 1— Explorer XII Energetic Particle Satellite. Figure 2— Geometry for solar radiation torques produced 

by a solar paddle. 
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( 3 ) 


Aj = A Q | cos 9 sin nZj sin <fi + cos nZj cos 4> | , 

A 2 = | “ sin 0 sin nz 2 sin <fi + cos nz 2 cos <f> | , 

A 3 - A 0 | - cos 0 sin nZj sin + cos nz x cos <j> | , 

A 4 = A 0 | sin 9 sin nz 2 sin 4> + cos nz 2 cos 4> | , 


where nz p nz 2 are the angles from the normal to the paddle to the spin axis and A 0 is the geometric 
area of the paddle. Absolute signs are used since a minus area means that the back side of a paddle 
is illuminated. Substituting Equation 3 into 2, Equation 2 into 1, and making the abbreviation 


M - sin nz sin 4> , 
N - cos nz cos <£> , 


the equation of motion becomes 

= PA 0 K(*,6>){- | 

+ | 
+ 1 


Mj | cos 9 | + Nj | hj | cos 9 | ^ 

- M 2 I sin <9 | + N 2 j h 2 | sin 9 

- Mj | cos 0 | + Nj | hj | cos 9 
M 2 | sin 6 | + N 2 J h 2 | sin 9 \ ^ 


sin + p sin(<£ - 2 nSj)^ 

| ^sin $ + p sin (<£ - 2 ns 2 )j 
| ^sin 4> + p sin((^) - 2ns 3 )^ 
sin 4> + p sin(^ ~ 2ns 4 


(4) 


Again absolute signs are used inside the parenthesis because there is 180 degree symmetry with the 
two paddles. That is one paddle presents the same area as the other paddle 180 degrees later. A 
shadow factor K is included to account for paddle on paddle shadowing and shadowing by the 
satellite structure. 

Strictly speaMng Equation 4 should contain a torque term due to albedo radiation from the earth. 
Because of the complicated mathematical description and the fact that it contributes very little in 
this application (it decreases as 1/a 2 where a is the ratio of the satellite distance to earth radius, 
Reference 3), the albedo radiation term will be neglected. 

As was noted earlier, the solar radiation factor, p, in Equation 4 is variable; it varies inversely 
as the square of the earth-sun distance. By considering the earth’s orbit around the sun as a per- 
turbed circular orbit, it can be shown that the distance between the sun and the earth as a function of 
time can be expressed as 


r(t) 


r o 


+ e cos 277 



( 5 ) 


where 

e = eccentricity of the earth’s orbit = .01674, 
T = 365 days, 

r 0 = mean earth- sun distances. 
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Time zero is at the maximum earth-sun distance of approximately June 21st. P can then be ex- 
pressed as 

P(t) ~ (1 + e cos 2ut/T) 2 

where p 0 is the mean solar constant at the mean solar-earth distance. P(t) is considered zero when 
the satellite is eclipsed by the earth. 

The change of satellite spin is obtained by integrating Equation 4; this is indicated symbolically 
in two ways below: 


CO - CO 


0 


[ P(t) 
J o 


}dt 


(7) 


or 


2A_ C 6 

- r P(t) K(4>,9) { }d0 , 


( 8 ) 


where the ns angles are obtained from the following equations 


COS 

ns j 

1 M, | 

cos 6 

1 + 

Nj 

cos 

ns 2 = 

1 - m 2 

| sin 

& 1 

+ N 2 

cos 

ns 3 = 

1 - M, 

| cos 

9 1 

+ N, 

cos 

ns . - 

4 

1 m 2 | 

sin 0 

1 + 

n 2 1 


(8a) 


and the quantity inside the brace is the same as that given in Equation 4. To carry out the integra- 
tion indicated by Equation 7 or 8 is impossible by analytical means and very messy and tedious to do 
numerically; a digital computer solution is required. This has been done and the results are dis- 
cussed below. It is of value though to develop an approximate equation (and solution) so that "ball 
park" answers can be determined readily. Within certain limits, it has been found that the approx- 
imate solution gives answers comparable to the digital solution. 


APPROXIMATE METHODS 

An approximate expression will now be developed for the torque equation. Directing attention to 
Equation 4 and noting that if ns, ns 2 , ns 3 and ns 4 are replaced by an average angle, <ns> , Equation 4 
can be put into the following form: 
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I dt " PA 0 K(0, ^) | ( | Mj | cos 0 | + N, | “ | Mj | cos 9 | + Nj tq | cos 6 | 

+ ( | - M 2 | sin 9 j + N 2 | — j M 2 | sin 0 | + N 2 I) h 2 ) sin 0 || (sin 4> + p sin (<£- 2 <ns> )) . (9) 

For 9 near zero and M > N (this restricts <p around 90°), the quantity in the brace in Equation 9 has 
the approximate value 

“ 2Njhj | cos 9 | 


and for 6 near 77/2 the brace has the approximate value 

- 2N 2 h 2 | sin 0 | 

For an approximate solution, it does not seem unreasonable to replace the complicated expression 
inside the brace by the sum of these two quantities 


“ 21^^ | cos 6 | - 2N 2 h 2 | sin 6 | 

throughout the paddle array rotation. Therefore, the following approximation can be made to Equa- 
tion 1. 


A° r 

-2-r j p(t)K<0,0)| 


cos nz 1 hj 


+ cos nz 2 h 2 




(sin + p sin (<4 ~ 2 <Cns> )jcos 4> dt . (10) 

The same result can be obtained by squaring a typical quantity | Mcos u + N | , neglecting n 2 , then 
taking the square root, expand in a power series and retain the first two terms. 

The average angle <ns> must now be evaluated. In an approximate equation it is desirable to 
have the equation as simple as possible and still retain the essence of the phenomena it is to describe. 
An expression can be developed for the average angle but it is complicated and unwieldy. The sim- 
plest value for <ns> would be zero. At first glance, this may appear to be an unwarranted simpli- 
fication. Subsequent comparison of numerical results from the theoretically correct equation will 
show that this is justified. Essentially, when <ns> is given a value of zero, it is tacitly assumed 
that the reflected light reinforces the direct light effect. 


Using the method of Kryloff and Bogoliuboff (Reference 4), that is, averaging over a cycle while 
considering du/dt constant during the averaging period, Equation 10 can be simplified to 


V - co 0 



R(t) K (cp) sin dt , 


(ii) 
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where 


R(t) 


1 

(1 + e cos 277 1 /T) 2 


Vo / 

1 + / : ) 77 1 V 1 ! 


h 2 cos nz 2 ! 


Equation 11 still cannot be integrated analytically since K(g*j) and <£ are usually known only graphi- 
cally. Using the trapezoidal rule, the integration can be replaced by a finite summation 

^ -“o = -PL S . R ( l .) K (^) sin 2 0 At ' (12) 

; = n 


where S fl = s N - 1/2, and the other S’s are equal to 1. K <p is known day by day, the increment of 
integration, At, can be a day without any appreciable loss of accuracy. For a quick estimate of 
solar radiation effect on the spin, K (</>) can be set equal to unity. 


DISCUSSION 

The exact spin equation (Equation 8) was applied to the Explorer XII satellite with the following 
parameters: 


A Q = 1.75 ft 2 , I - 3.53 slug ft 2 , 

h, = h 2 " 2.35 ft, P 0 = 9.50 x 10~ 8 lb/ft 2 , 

nz j “ nz 2 = 59° 47*, P w 0.10. 

“i = M 2 . n i = N 2> 

The reflectivity coefficient for a solar paddle is a rather uncertain quantity; the above value for 
was estimated from tests conducted on plain solar cells. A plot of the telemetry optical aspect 
sunline-spin axis angle used in the calculations is shown in Figure 3. The shadowing factor, K vs #■ , 
is shown in Figure 4. Since Explorer XII was in the sunlight more than 98 percent of the time, no 
eclipse factor was used. Equation 8 was solved on a 7090 digital computer in an iterative procedure 
by calculating the change of spin over one rotation and considering this change of spin applied for 
200 revolutions. The elapsed time was calculated and the sunline-spin axis angle was changed ac- 
cordingly; the calculation was then repeated. The above method resulted in reasonable computer 
times and still gave, it was felt, accurate results. Computer results are compared with actual sat- 
tellite spin data in Figure 5. Calculations were started two days after launch to allow the nutation 
angle to dampen out since after yo-yo operation (a spin reduction device), Explorer XII had a five 
degree half cone angle nutation. The optical aspect sensor indicated that after the initial coning 
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SHADOWING FACTOR 








Figure 5— Spin rate of Explorer XII vs time. 


angle damped out, there was virtually no nutation angle during the active life of the satellite. It will 
be noticed from the figure that there is good agreement between the theoretical curve and actual spin 
data. Maximum deviation from the theoretical curve is about one percent. Considering the uncer- 
tainty of the reflectivity coefficient and the small but finite perturbations of the albedo, magnetic, 

and aerodynamic torques, the agreement is satisfactory. 

The approximate equation given by Equation 12 was also applied with the Explorer XII para- 
meters (At was taken to be one day) and results are compared with the computer solution. This is 
shown in Figure 6. It will be noted that there is very little difference between the computer solution 
and the approximate solution. For most practical purposes the curves are essentially the same. Be- 
cause of the good agreement it is inferred that the approximate equation gives satisfactory results 
for 0 between 45° and 135° and cos nz < 1/2. 

It is of value to note that in general solar radiation torques about the spin axis occur because of 
the different tilt (projected area) of the paddles. This is clearly seen in Equation 2; if the projected 
areas are equal, (ns angles are also equal) there are no solar torques about the spin axis. For small 
reflectivities, it can also be stated that if the sun shines from below the equator of the satellite 
(<£ > 90°), the satellite will spin up; if the sun shines from above the equator (<£ < 90°), the satellite 
will spin down. This is easily seen from the approximate Equations 10 or 11. This effect has been 
observed in satellites launched after Explorer XII. It is interesting to note that the maximum value 
of the solar torque around the spin axis is about 2 * 10 ‘ 7 ft. lbs. This is a very small torque, indeed. 
As was noted though, a small but persistent torque acting on a satellite for a long period of time does 
have an appreciable effect. 
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Figure 6— Explorer XII spin rate vs time: comparison of computer calculation with 
approximate solution (Equation 12). 


RESUME 

An equation of motion about a satellite spin axis has been developed that includes spectral solar 
radiation torques. The equation was applied to the Explorer XII to account for the observed increase 
in spin. Calculated spin values agree within about one percent of the observed satellite spin values. 
It is concluded from the analysis and calculations that Explorer XII spin-up was due to solar radia- 
tion torques. Also, an approximate spin equation was developed which permits rapid calculations of 
the effect of solar radiation on satellite spin. 

ACKNOWLEDGMENT 

The author is indebted to W. R. Mentzer of GSFC for performing the hand calculations, drawing 
the curves and figures and programming the equations on the IBM 7090 Computer. 


REFERENCES 

1. McElvain, R. J., "Effects of Solar Radiation Pressure upon Satellite Attitude Control,” in: Pro- 
gress in Astronautics and Rocketry edited by Roberson, R. E., and Farrior, J. S., New York: 
Academic Press, 1962. 

2. Clancy, T. F., and Mitchell, T. P., "Effects of Radiation Forces upon the Attitude of an Artifi- 
cial Earth Satellite, " Center for Radiophysics and Space Research Report #129, September 1962, 
Cornell University, Ithaca, New York. 




3. Levin, E., "Reflected Radiation Received by an Earth Satellite", ARS Journal 32(9) : 1328- 1331, 
September 1962. 

4. Kryloff, N., and Bogoliuboff, N., ’’Introduction to Non-Linear Mechanics” (translation), Princeton. 
Princeton University Press, 1947, chapters 10-12. 

5. Curtiss, D. R., and Moulton, E. J., "Essentials of Trigonometry", New York: D. C. Heath and 
Company, 1942. 


11 



Appendix A 


Projected Areas of a Solar Paddle 


The effective projected area of a solar paddle will now be determined. We wish to derive an ex- 
pression involving the angle between the spin axis of the satellite and the sunline, 0 , the rotation 
angle 8 , and geometrical properties of the solar paddles. Initially the paddle will be considered 
stationary and any shadowing by the satellite will be neglected. 

Let n be the unit normal to the paddle and 7 the unit vector in the direction of the sun (Figure Al). 
The projected area perpendicular to 7 is 


A 1 A n cos ns 
p 0 


(Al) 


where 

A = effective paddle area, 

p 

A 0 = geometric paddle area, 
ns = angle between vector n and 7. 


Now cos ns can be obtained from the vector dot product 


cos ns - n • s . 


(A2) 


n 



PADDLE 


If we choose the coordinate system shown in Figure A2, where 
i, j, k are unit vectors, the n and s vectors have the following 
components: 

n components - cos nx, cos ny, cos nz, 

7 components - sin 0 , 0, cos <f>. 

Carrying out the dot product results in 

n * s - cos nx sin 0 + cos nz cos 0 . (A3) 

Hence, Equation Al can be written as 

A p = A fl (cos nx sin 4 > + cos nz cos $) • (A4) 
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Equation A4 is for the static case (no spin). When spin is 
present we have the situation shown in Figure A3. Though sis 
not shown in the figure for clarity, it is still pointing in the di- 
rection of the sun. The z axis is considered to be the spin axis 
and therefore the angle nx will vary with time (see Equation A4). 
From solid trigonometry we can write the following formula 
(Reference 5, pp. 112-113) 

cos nx = cos 6 cos(rr/2 " nz) - cos 0 sin nz . (A5) 

Substituting Equation A5 into Equation A4 results in 


Figure A3 


sin nz sin 4 > cos 0 + cos nz cos 4 > 


Absolute signs are used since a negative area means that the back side of a paddle is illuminated. 
Equation A6 is the basic equation for the projected area of a solar paddle. 
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